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Abstract 



O We have developed a new technique to measure vis- 
j » coelasticity in soft materials such as polymer solutions, 
& by monitoring thermal fluctuations of embedded probe 
£3 particles using laser interferometry in a microscope. In- 
terferometry allows us to obtain power spectra of fluc- 
Ch tuating beads from 0.1 Hz to 20 kHz, and with sub- 
^ nanometer spatial resolution. Using linear response the- 

i iory, we determined the frequency-dependent loss and 

storage shear moduli up to frequencies on the order of 
a kHz. Our technique measures local values of the vis- 
1—1 coelastic response, without actively straining the sys- 
^ tern, and is especially suited to soft biopolymer net- 
Q\ works. We studied semiflexible F-actin solutions and, 
as a control, flexible polyacrylamide (PAAm) gels, the 
latter close to their gelation threshold. With small parti- 
cles, we could probe the transition from macroscopic vis- 
^ coelasticity to more complex microscopic dynamics. In 
— the macroscopic limit we find shear moduli at 0.1 Hz of 
^ 6" = 0.11±0.03 Pa and 0.17±0.07 Pa for 1 and 2 mg/ml 
Ch actin solutions, close to the onset of the elastic plateau, 
q knd scaling behavior consistent with G*(uj) ~ cj 3 / 4 at 
higher frequencies. For polyacrylamide we measured 
plateau moduli of 2.0, 24, 100 and 280 Pa for crosslinked 
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gels of 2, 2.5, 3 and 5% concentration (weight/volume) 



d respectively, in agreement to within a factor of two with 
values obtained from conventional rheology. We also 
found evidence for scaling of G*(oj) ~ w 1 / 2 , consistent 
with the predictions of the Rouse model for flexible poly- 
mers. 

Introduction 

The usefulness of synthetic polymeric materials, and the 
functions of biopolymers, are all largely based on com- 
plex mechanical properties that derive from a hierarchi- 



cal structure. The mechanical response of polymer solu- 
tions or gels displays characteristics of both fluids (vis- 
cosity) and of solids (elasticity), depending on the rate 
of change of applied stress. This response is convention- 
ally described by frequency-dependent storage and loss 
shear moduli which are commonly measured by active, 
mechanically-imposed oscillatory strain in macroscopic 
samples, in contrast to stationary-flow geometries used 
for measurements in fluids 1-5 . 

Such methods have also been used in the past to study 
reconstituted biopolymer networks in vitro. A gen- 
eral distinction between synthetic polymers and biopoly- 
mers is that the former are typically flexible, since their 
monomer size is small, whereas many biopolymers are 
formed from large protein monomers and thus are much 
less flexible. The viscoelastic properties of flexible and 
semiflexible polymer systems are quite different. At a 
given volume fraction of polymer, the shear modulus (or 
stiffness) of a semiflexible polymer network can be sev- 
eral orders of magnitude larger than that of a flexible 
one. This may be one reason why biological evolution 
has favored semiflexible polymer networks for mechani- 
cal stability. F-actin is one of the primary components of 
the cytoskeleton of plant and animal cells, and is largely 
responsible for the viscoelastic response of cells 6,7 . The 
~17 //m persistence length of actin filaments 8,9 is about 
three orders of magnitude larger than their diameter. 
This tremendous aspect ratio makes them ideal model 
semiflexible polymers. They are rather rigid on the scale 
of cytoskeletal networks of typical cells, which have char- 
acteristic mesh sizes on the order of microns. In recon- 
stituted F-actin networks, individual filaments can be 
hundreds of microns long, so that solutions of less than 
0.1% polymer (volume fraction) are still strongly entan- 
gled. 

Macroscopic rheological measurements in reconsti- 
tuted F-actin solutions have been reported in the 
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literature . These experiments have shown that the 
regime of linear viscoelastic response is small, 16 of or- 
der a few percent, and that this range decreases with in- 
creasing F-actin concentration. For larger strains, strain 
hardening is observed, followed by apparent shear thin- 
ning, which may be due to filament breakage. This is 
to be expected for semiflexible polymer networks, since 
the amplitudes of thermal undulations are small, and 
the full extension of individual filaments is thus reached 
for small strains, 20 making the response nonlinear. This 
makes measurements with commercial rheometers prob- 
lematic, although some custom intruments can deal with 
strains substantially smaller than Active macro- 

scopic measurements on actin have given controversial 
results 16 and nonlinearity may be one reason for this. 

In attempts to approach in vivo conditions, viscoelas- 
ticity has also been studied on a microscopic, cellular 
scale. Early work included the manipulation of mag- 
netic particles in gelatin 21 and in living cells 22 . Recent 
experimental techniques for manipulating and tracking 
sub-micrometer particles have revived the interest in mi- 
crorheology. Several studies have used direct manipu- 
lation of embedded magnetic beads in solutions of F- 
actin 23-27 . These experiments have been limited by the 
spatial and temporal resolution of video microscopy. 




Figure 1: Schematic of a probe particle of radius R em- 
bedded in an actin network of mesh size £. 

We introduce here a microscopic method of measur- 
ing viscoelastic properties within micrometer-sized sam- 
ple volumes. We study entangled, but non-crosslinked 
networks of semiflexible F-actin, as well as low vol- 



ume fraction gels of crosslinked polyacrylamide (PAAm) , 
a flexible polymer. We observe the thermal fluctua- 
tions of micrometer-sized particles embedded in soft gels 
(Fig. 1). Laser interferometry in a light microscope pro- 
vides high resolution (less than 1 nm) and bandwidth 
(from 0.1 Hz to 20 kHz). Using dispersion relations from 
linear response theory, the frequency-dependent loss and 
storage shear moduli can be determined from the fluc- 
tuation power spectra. The technique is in principle less 
invasive than active methods, in that no strain at all is 
imposed on the material. 

Diffusing wave spectroscopy (DWS) has been used by 
others to observe thermal fluctuations of ensembles of 
particles, and viscoelastic properties in polymers and 
colloids have thereby been deduced 28,29 . This method 
measures average viscoelastic properties, as opposed to 
our local measurements; we further compare this con- 
ceptually related method with ours in the Discussion. 

Our method measures shear moduli over a larger fre- 
quency range than accessible to video based microrheol- 
ogy (while DWS and some other macroscopic methods 
can reach higher frequencies 1 ' 30 ). Besides providing a 
new way of passively measuring shear moduli in both 
synthetic and biological polymer systems, our technique 
can be sensitive to dynamics and material parameters 
that are inaccessible to macroscopic methods: spatial 
inhomogeneities in networks can be studied on a mi- 
crometer scale. At low frequencies, our method appears 
to be sensitive to dynamics, such as non-shear deforma- 
tions of the network with respect to the solvent, that 
are not seen in macroscopic mechanical rheology. For 
probe sizes comparable to the mesh size of the polymer 
network, deviations from continuum elasticity become 
apparent. With decreasing probe size, the transition 
from collective network dynamics to single-filament dy- 
namics can be traced to explore the microscopic basis 
for macroscopic properties. This is particularly relevant 
for biopolymer networks, with mesh sizes as large as mi- 
crons. More importantly for biology, microrheology will 
permit the characterization of small samples, such as 
living cells. 

Theory and Data Analysis Methods 31 

We model our experiments as embedding a spherical 
bead, of radius R, in a linear viscoelastic medium of 
density p with a macroscopic complex shear modulus 
G* and Poisson ratio v? 2 For motions of micron-sized 
probe particles, inertial effects are negligible at frequen- 
cies uj where the inertial decay length (i.e., the inverse 
magnitude l/|r| of the complex propagation constant 1 
r = (pu? /2G*) 1/2 ) is lar ger in magnitude than R. At 
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1 kHz, 1/|T| is at least tens of microns in all our sam- 
ples. Correspondingly, inertial effects can be expected 
to be relevant only for frequencies of order 1 MHz. 
In addition, the Reynolds number for thermal motions 
on this scale is small: ~ 10~ 3 . In an incompressible 
medium, bead motion is then determined entirely by the 
frequency-dependent complex shear modulus G*(uj) = 
G'(uj) + iG"(u), where G' and G" are the storage and 
loss moduli. 

In extremely low-frequency motions of a gel, the net- 
work deforms so slowly that stress in the solvent relaxes 
completely, and the network can undergo not only shear 
but compressional deformation, i.e. it behaves as an 
elastic continuum with shear modulus G = G*{u) = 0) 
and Poisson ratio v < 1/2. Stress cr^ is related to strain 
Uij by 32 



aij = 2G [uij + i/<5jjEu fcfe /(l-2i/)] 



(1) 



The equation of elastic equilibrium, d<Jij/dxi = 0, can 
be solved exactly for with no-slip boundary condi- 
tions on a rigid spherical surface, 33-35 yielding an effec- 
tive compliance (reciprocal spring constant) for sphere 
displacement 
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The static compliance a in the limit u — > is significant 
because it is related by equipartition of energy to the 
total mean square fluctuations of the probe particle in 
the gel: 

(x 2 ) = kTa(io = 0). (3) 

Eqs. (3) and (2) together provide an estimate of the 
static modulus G from the variance (x 2 ) of fluctuations: 
a priori, v is unknown in Eq. (2) but over the relevant 
range, < v < 1/2, v does not strongly affect the re- 
sult; one can put a ~ \/£>itGR with at most 25% error. 
However, (x 2 ) is often dominated by low- frequency in- 
strumental drifts. We will instead use a method (to be 
described below) of estimating G*(uj = 0) that is less 
sensitive to such noise. 

For frequencies above the u — ► limit, viscous stresses 
develop in the gel, and we must reconsider Eq. (2). The 
solvent and the polymer network are coupled through 
viscous drag, which becomes stronger with increasing 
frequency. Eventually, solvent and network are expected 
to move as one at scales large compared with the network 
mesh size £. In this case the polymer network is strongly 
coupled to the incompressible solvent and so behaves like 
an incompressible network, which can be described by a 
Poisson ratio v = 1/2. The crossover frequency, above 



which incompressible behavior is observed, can be esti- 
mated as follows. The viscous force per volume exerted 
by the solvent on the network 36 ' 37 is ~ f]v/^ 2 , where v is 
the velocity of the solvent relative to the network. The 
local elastic force per volume (exerted by the rest of the 
network) is GV 2 u ~ u/R 2 at the bead surface, where u 
is the network displacement field. Viscous coupling will 
dominate above a crossover frequency 
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Eq. (4) is an order-of-magnitude estimate of ui c that 
corresponds to a few Hertz in our F-actin solutions 
[G ~ 1 Pa, £ ~ 0.1R, rj = 1 cP). Above co c , provided 
that £ is small compared with R, the bead moves as 
an inclusion in an incompressible continuum viscoelas- 
tic medium. The stress-strain relations in Eq. (1) still 
hold, but with v = 1/2. The spatial equation of motion 
is still daij/dxi = (at a given id) and its solution is 
Eq. (2) with v = 1/2: 
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a*(u) and G*(u>) are complex because this formula de- 
scribes sinusoidal sphere displacements x(t) = a*{uj)f{t) 
caused by a force f(t) oc e -lw *. This is a generaliza- 
tion of the well-known Stokes formula / = Qirqrx, since 
rj = iG*(uj)/uj in a purely viscous fluid. 

We can rigorously justify Eq. (5) directly, as follows. 
It applies to an incompressible medium without inertia. 
The dynamics of the medium are completely determined 
by its stress-strain relation (Eq. (1)), which becomes 
&ij = 2G*{ui)uij. This is mathematically the same law 
that describes a simple viscous fluid, but with the (com- 
plex and frequency-dependent) quantity iG*{uj)/oj re- 
placing the viscosity rj. Together with the equation of 
motion, daij/dxi = 0, the spatial equations to be solved 
are also mathematically identical. The unique solution 
to these equations at a frequency uj, that satisfies no- 
slip boundary conditions on the surface of the moving 
sphere, is precisely the Stokes law, with the (complex 
and frequency-dependent) G*{ui) replacing the quantity 
—iu>r). This has the simple and remarkable consequence 
that the well-known Stokes force generalizes to give the 
correct compliance, involving the complex G*{u). 

For u between zero and u c , the solvent and network 
may undergo significant relative motion. As a result 
of this partial "draining" of the network, Eq. (5) will 
not be valid; the situation must then be described by a 
"two-fluid" model 36,3 . Force response and fluctuations 
may become dominated by long-wavelength modes that 
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relax over distances large compared with ^(G/ijuj) 1 ^ 2 . 
This low-frequency draining allows for relative motion of 
network and solvent, involving not only shear modes but 
also compressional modes in the network. As the spatial 
extent of the modes becomes larger than the bead size, 
the bead motion and its power spectrum are expected 
to become independent of bead size. 

The compliance a*(io) = a'(u) + ia" '(u) is the com- 
plex response function for bead displacement. Thus, 
for frequencies above id c , Eq. (5) and the fluctuation- 
dissipation theorem 38 relate (x 2 ), the power spectral 
density (PSD) of thermal motion, to the imaginary part 
a" {id): 

2 _ 4kTa"(io) 
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Provided that a" (to) is known over a large enough fre- 
quency range, one can recover the real compliance a'{oj) 
from the Kramers-Kronig relations 38 by evaluating a dis- 
persion integral 
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This allows us to explicitly calculate (by Eqs. (5) and 
(6)) the storage and loss moduli G'{u) and G"(lj). P 
in Eq. (7) denotes a principal-value integral, meaning 
the e — > limit of the sum of two ^-integrals: from 
to {to — e), and from {u + e) to oo. However, as writ- 
ten in Eq. (8), this is equivalent to successive sine and 
cosine transforms of a"(co), which are conveniently per- 
formed for the discrete data points of a long time series. 
Standard Fourier routines 39 can be used to evaluate the 
discrete PSD and then to perform the sine and cosine 
transforms. Either the PSD or the output of Eq. (8) can 
be smoothed (by averaging within bins of equal logarith- 
mic spacing) without strongly affecting the numerical 
results. As a control, we have also calculated Eq. (7) di- 
rectly in the logarithmic domain, finding similar results. 

When the behavior of the network plus solvent is well 
described as a single-component medium (so that Eq. (5) 
holds), we can take the complex reciprocal of a* (to) to 
obtain the complex G*(id): 
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Qualitatively, we expect the shear modulus G*(oj) to 
exhibit a characteristic form, 4 involving up to three dis- 
tinct dynamical regimes. If the polymer network is en- 



tangled but not crosslinked, as in F-actin solutions, vis- 
cous flow will occur over times longer than the repta- 
tion time r r , with G'(u) going to zero. The reptation 
time is hours or days, for actin filaments tens of microns 
long, 40 and is thus not accessible to our experiments. 
For frequencies above l/r r , a rubber- like plateau ap- 
pears, with a frequency-independent elastic response G', 
as in a crosslinked gel. Above the high-frequency end of 
the plateau, the moduli are expected to increase with a 
characteristic power of frequency, G' and G" ~ oj z , re- 
flecting the increasingly limited relaxations of dynamic 
modes within a mesh of the network 4 . Power law be- 
havior is expected for the shear moduli of any polymer 
network above the characteristic mesh relaxation time 
and below the molecular high frequency cut-off, since 
there is no other characteristic length or time scale in 
this regime. In this scaling regime, from Eqs. (5) and 
(6), the PSD is then expected to follow 
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(7) Materials and Experimental Methods 



Actin was purified (^95% purity) from chicken skele- 
tal muscle, following standard recipes 41 . Its con- 
centration was determined both by staining (BioRad) 
and by UV absorption at 290 nm (specific absorption: 
0.65 cm 2 /mg). Monomeric actin (G-actin) was rapidly 
frozen in liquid nitrogen and stored at — 85° C. Samples 
were prepared by mixing G-actin with a small number of 
silica beads (Bangs Laboratories, except for the 0.5 /jm 
diameter beads, which were kindly provided by E. Mati- 
jevic). After adding concentrated polymerization buffer 
(F-buffer: 2 mM hydroxyethyl-piperazineethanesulfonic 
acid (HEPES) (pH 7.2), 2 mM MgCl 2 , 50 mM KC1, 
1 mM ethylenebis(oxyethylenenitrilo)tetraacetic acid 
(EGTA), 1 mM adenosinetriphosphate (ATP)) the mix- 
ture was immediately transferred into a sample chamber 
made from a microscope slide, a cover slip and double- 
stick tape (with inner dimensions 15mm x 3mm x70//m). 
In the sealed sample chamber the actin polymerized at 
room temperature for at least one hour under slow ro- 
tation. The samples were stored at 4°C and examined 
within one day (and in one case remeasured the following 
day, as a control). 

Polyacrylamide (PAAm) gels were prepared accord- 
ing to a standard gel electrophoresis recipe, 42 with 
concentrations of 2, 2.5, 3 and 5% (weight/volume) 
and a relative concentration of 3% bis-acrylamide as 
a crosslinker. Solutions were thoroughly degassed un- 
der vacuum, beads added and polymerization initiated 
with tetramethylethylenediamine (TEMED) and ammo- 
nium persulfate (APS). As with actin, polymerizing 
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solutions were then transferred into sample chambers 
(15mm x 6mm x 140/um) and slowly rotated at room tem- 
perature for at least one hour before starting experi- 
ments. 
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Figure 2: Schematic diagram of the custom-built in- 
verted light microscope used for differential interference 
contrast imaging and laser interferometric displacement 
detection. 

Our microscope is a custom-built inverted instrument 
(optics from Carl Zeiss, Inc.), constructed on an optical 
rail system and mounted on a vibration isolated optical 
bench (Fig. 2). To detect the thermal motion of beads 
imbedded in the gel we used an interferometer 3 ' (see 
Fig. 3) with near- infrared laser illumination (1064 nm 
NdVC-4, 3.4 W (cw) max. power, Topaz 106c, Spectra 
Physics) . 

A linearly polarized laser beam is split in two beams 
by the Wollaston prism below the objective, which pro- 
duces two diffraction limited overlapping foci in the spec- 
imen plane. A refractive particle located asymmetrically 
within the two foci will cause slightly elliptical polariza- 
tion after recombination of the two beams by the upper 
Wollaston prism. A quarter- wave plate renders the light 
close to circularly polarized, whose two perpendicular 
linear components are then detected by two photodi- 
odes. Deviations from circularity are calibrated to mea- 
sure particle displacements. The normalized difference 
between the two signals is calculated by custom-built 
analog electronics (Fig. 3), amplified and anti-alias fil- 
tered slightly below the Nyquist frequency. This analog 
signal is then digitized and recorded using an A/D in- 
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Figure 3: Schematic diagram of the interferometric dis- 
placement detection system. 



terface (MIO 16X, National Instruments) and data ac- 
quisition software written in Labview (National Instru- 
ments) . 

Although this is not done here, the focused laser beam 
can act, at high enough laser power, as an optical trap, 
exerting forces on the particle in the focus 44,45 . Here, 
we want to measure unperturbed thermal fluctuations, 
and thus need to minimize the trapping force. For that 
purpose, a 1.5 mm diameter pinhole was inserted at a po- 
sition conjugate to the back focal plane of the objective, 
broadening the laser focus and increasing the detector 
range. The laser power was typically 0.6 mW in the 
specimen, low enough to make optical forces negligible. 

Before recording fluctuations of an individual bead, 
it was centered in the detector range using a piezo- 
actuated translation stage. For all experiments reported 
here, the data acquisition rates were 50 or 60 kHz (anti- 
alias filtered at 25 kHz). Time series data were recorded 
for at least 17.5 seconds (2 20 data points) to obtain 
power spectra ranging from 0.1 Hz to 25 kHz. Each spec- 
trum was evaluated by Fast Fourier Transform (FFT) 
after applying a Welch window. Power spectra were 
smoothed by averaging within bins of geometrically in- 
creasing width (with factor 1.1). 

The linear range of the detector is about 200 nm for 
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Figure 4: Detector sensitivity for displacement of silica 
beads of various diameters embedded (A) in an actin so- 
lution (2 mg/ml F-actin) and (B) in a 2.5% PAAm gel, 
determined by driving a piezo-actuated stage with a tri- 
angular signal at a frequency of 0.1 Hz. Response curves 
are averages over several periods, sampled at 100 Hz (in- 
set), smoothed (by a Savitsky-Golay filter) and shifted 
vertically to pass through the origin. For a uniform scale 
some curves include a gain factor. Inset: Direct signal 
time series for a 0.5 /xm bead. 
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Figure 5: Bead-size dependence, and scatter, of detector 
response factors (see Fig. 4) . Four different experiments 
are compared; two with F-actin and two with PAAm. 
Data points are horizontally shifted for clarity. Vertical 
scatter is partly due to bead polydispersity, but PAAm 
experiments show a distinctly stronger scatter. Bead 
diameters are 0.5, 0.9, 2.1, and 5.0 /urn. 



0.5 jira beads and increases with bead size. Each ob- 
served bead was immediately calibrated following the 
recording of thermal motion, by moving the sample on 
the piezo-actuated stage with a triangular signal wave 
form, producing constant velocity displacements (be- 
tween 150 nm and 1 /j,m peak-to-peak) at a frequency of 
0.1 Hz. The driving voltage was digitally synthesized to 
correct for the non-linearity and hysteresis of the piezo 
actuators. Driver voltage (not shown) and detector re- 
sponse (Fig. 4 inset) were recorded for about 40 seconds 
and analyzed by plotting detector response versus bead 
displacement. Slopes were estimated by linear regression 
from an averaged response-displacement curve for each 
individual bead (Fig. 4A and B). The result is a sen- 
sitivity factor which is used to convert detector signals 
into actual displacements. 

While centering the bead in the detector range, we 
found that the sensitivity does not vary strongly in the 
perpendicular direction (Y), out to Y-offsets compara- 
ble to the linear range. The sensitivity also did not vary 
significantly with axial position in a range for which the 
bead appeared focused in the video image (~ 0.6 fim). 
Detector sensitivities as a function of bead diameter are 
plotted in Fig. 5. They increase roughly with the third 
power of bead radius as long as the bead is smaller 
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than the laser focus, since the scattering amplitude is 
expected to scale with the volume of the scatterer. For 
larger beads (in the limit of geometrical optics), the sen- 
sitivity decreases roughly as i? -1 , the only remaining 
length scale. In other words, the same relative displace- 
ment Ax/R should give the same signal AV, indepen- 
dent of bead size: AV oc Ax/R. 

Investigating the origin of scatter in the detector sensi- 
tivity, we found that variations in beam alignment dom- 
inated bead polydispersity for the actin experiments. 
Bead polydispersity was checked directly with transmis- 
sion electron microscopy (TEM; data not shown) and 
found to be about 10%, with a non-normal distribu- 
tion. The images showed both larger beads, which ap- 
pear to have formed by merging nuclei, and very small 
beads that may be contamination. The scatter in the 
displacement-response slopes was about ±10% (coeffi- 
cient of variation) for actin. When no correlation was 
evident between calibration factors and the amplitudes 
of the PSD at around 100 Hz (which should scale with 
bead radius R) we averaged the factor over all beads 
observed in a sample. 

For PAAm the scatter in sensitivities for a given bead 
size appeared significantly larger (±30-60% coefficient of 
variation), and response factors were indeed correlated 
with the PSD, indicating that at least part of the scatter 
was caused by bead size polydispersity. In these cases, 
spectra were calibrated by the factors determined for 
individual beads, which decreased the variance in the 
spectra between beads in the same sample. Since the 
response factor scales roughly with the third power of 
bead diameter for the smaller beads (see Fig. 5), a 10% 
coefficient of variance can explain a 30-40% scatter in 
the spectra. The polymer network itself and possible 
surface adsorption is not expected to influence apparent 
bead size and thus sensitivity substantially. Further- 
more, the index of refraction changes little with polymer 
volume fraction in PAAm (measured to be about a 1.2% 
increase in relative index of refraction at 10% polymer 
w/v). Finally, surface chemistry that could produce a 
dense layer of polymer directly on the bead surface is 
not expected to produce more than a monolayer, and is 
thus limited to a thickness on the order of tens of A. 

Results 
Actin solutions 

Fluctuation time series were recorded for silica beads of 
0.5, 0.9, 2.1, and 5.0 /jm diameter, embedded in actin 
gels of 1 and 2 mg/ml concentrations. The fluctuation 
signal is well above the noise, as shown by comparison 
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Figure 6: PSDs (power spectral densities) of background 
noise compared to signal for F-actin and PAAm. Volt- 
age signals were sampled at 50 kHz and spectra were 
smoothed by averaging in equal log bins. An anti- 
aliasing filter attenuates the signal strongly above about 
20 kHz. The lower two curves show the PSD of the 
interferometer signal for water without bead, and with 
a 0.5 /Ltm bead fixed on a glass surface. The F-actin 
spectrum (0.5 //m bead in 2 mg/ml actin) is everywhere 
about two orders of magnitude above the noise, whereas 
the PAAm spectrum (0.5 /im beads in 2% PAAm) is 
dominated by noise below about 1 Hz. 

to spectra of a bead immobilized on a surface, and of the 
water background without bead (Fig. 6). These controls 
display the effects of low-frequency instrumental noise 
due to thermal drifts in the microscope and the laser 
beam path, acoustic vibrations, and beam-pointing fluc- 
tuations. At high frequencies, noise depends on laser in- 
tensity since detection shot noise and preamplifier noise 
are dominant there. We emphasize that the amplitude 
of this noise is small, and remains about two orders of 
magnitude below the signal. Note that all spectra are 
(analog) anti-alias filtered at 25 kHz, causing the "tails" 
at about 20 kHz. 

The power spectra of bead fluctuations in actin solu- 
tions were reproducible in shape and amplitude between 
different beads in the same sample (Fig. 7A), as well as 
between samples (Fig. 7B). Observed coefficients of vari- 
ance (of the PSD at 100 Hz in Fig. 7A) within a sample 
(N=6) were 22% (0.5 /im), 24% (0.9 /j,m), 8.7% (2.1 fim) 
and 15% (5.0 /jm bead diameter). Fig. 7B shows two 
spectra from samples in different experiments but from 
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the same protein preparation, and one made with actin 
from a different laboratory (sample provided by Paul 
Janmey). We find that beads in different samples do 
not scatter more than different beads within one sam- 
ple. 

Slow changes have been reported to occur in actin so- 
lutions over many hours after polymerization, possibly 
due to the depletion of ATP in the buffer, or to the an- 
nealing of filaments. To address such effects, we remea- 
sured samples the day following their preparation, with- 
out seeing significant differences in the spectra (Fig. 7C). 
Observed beads were typically located at heights be- 
tween 10 and 50 /j,m above the sample chamber surface. 
Variations in the spectra were not correlated with dis- 
tance from the boundaries (data not shown). 

Power spectral densities (PSDs) of thermal motions of 



Figure 7: (A) Reproducibility within individual samples 
of F-actin (2 mg/ml). Six beads each with diameter 0.5 
and 5.0 /j,m. Height above the sample cell surface was 
between 10 and 50 /im. (B) Reproducibility between 
experiments and different protein preparations. Spectra 
from two experiments using the same actin preparation 
at a concentration of 2 mg/ml (#1: 1.8 jum, #2: 2.1 /im 
silica beads) and using actin from a different laboratory 
(P. Janmey) (#3: 1.8 /im silica beads). (C) Aging of F- 
actin samples. Spectra from four samples (0.5, 0.9, 2.1, 
5.0 fim beads, top to bottom, with six beads averaged 
per sample) measured within hours of polymerization 
and remeasured the following day. (D) Comparison of 1 
and 2 mg/ml F-actin solutions. Spectra from two sam- 
ples with 0.9 ^m beads and two samples with 5.0 /urn 
beads. Six beads are averaged in each plotted PSD. 
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0.9 and 5.0 fim silica beads in actin solutions at different 
concentrations are shown in Fig. 7D. The shape of 1 and 
2 mg/ml spectra is similar and the PSD amplitude of the 
less concentrated sample is smaller, as expected. 

A phase transition to a nematic phase has been re- 
ported for actin at concentrations of about 2 mg/ml 40 ' 46 . 
However, we see no qualitative difference in the ap- 
pearance or in the behavior of our actin gels at 1 and 
2 mg/ml. There were no signs of inhomogeneities un- 
der polarization microscopy. Furthermore, we do not 
see any greatly increased scatter in our measured values 
of G*, which would be expected in a highly anisotropic 
medium, given our uniaxial detection technique. 

Fig. 8A shows averaged and smoothed spectra ob- 
tained from 2 mg/ml actin solutions with beads ranging 
from 0.5 to 5.0 /urn in diameter. The variance in the 
smoothed spectra (values at 100 Hz) as a function of 
bead diameter is plotted in the inset. Part of the scat- 
ter presumably reflects local inhomogeneities in the net- 
works. Smaller beads, which probe more locally, should 
show more scatter. For 0.5 /jm beads, however, some 
variation is also due to drift out of the linear range of 
the detector, and to bead polydispersity, as discussed in 
Materials and Experimental Methods. Detector sensi- 
tivity depends on bead size, as we have shown in Fig. 5. 
The high frequency noise level follows this dependence, 
with the smallest (0.5 fim) and largest (5.0 /xm) beads 
showing the largest relative noise levels. 

The largest beads are about a factor of 10 to 20 
larger than the average mesh size (0.25 /um for actin 
at 2 mg/ml 47 ) and a continuum elastic model should 
apply. Approximate power law behavior, with a slope 
of about —1.75, is evident over about three decades in 
frequency for the 5 /xm beads, as indicated in Fig. 8. All 
spectra show a slight downturn at about 3 kHz, an ef- 
fect which we cannot yet explain. At the high frequency 
end, the data have not reached the noise bottom evident 
in the spectrum of a bead fixed to the substrate surface 
(Fig. 6). 

Eq. (5) predicts that spectra should scale with the 
bead radius R for high frequencies. This is roughly ob- 
served as shown in Fig. 8B, but a systematic deviation 
from i?-scaling is evident. Small beads show power spec- 
tral densities of higher amplitude than expected. We be- 
lieve that the increased fluctuations are not due to dif- 
fusion through the network or constraint release due to 
filament reptation; in video recordings of up to two hours 
we did not observe any long-range diffusive motion. Fur- 
thermore, we can exclude that the beads themselves bind 
to F-actin (with potential bunching of filaments as a re- 
sult), because we observed practically unimpeded diffu- 
sion of silica beads (0.2 /im) with diameters below the 
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Figure 8: (A) Bead-size dependence of PSDs in 2 mg/ml 
F-actin solutions using a 60 kHz sampling rate. Spectra 
were calibrated with response factors found as in Fig. 4. 
PSDs from beads between 0.5 and 5.0 /im diameter are 
shown. For each bead size PSDs from six different beads 
in one sample were averaged. The line indicates a power- 
law slope of —1.75. Inset: PSD amplitudes at 100 Hz 
of individual smoothed spectra that make up the aver- 
ages. (B) Averaged PSDs from (A) multiplied by their 
respective bead diameters. Inset: direct time series data 
for 0.9 /j,m bead. 
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mesh size. 

We instead attribute the lack of scaling to the steric 
exclusion of filaments near the bead. This model is 
elaborated quantitatively in a separate publication 34 but 
given qualitatively in the following. A large bead of ra- 
dius R will affect the distribution of filaments to an ex- 
clusion depth I ~ £ p , where l v is the filament persistence 
length. Since the Stokes flow field around a large bead 
extends beyond the exclusion depth to a distance on the 
order of R, large beads will mainly see the unperturbed 
network according to Eq. (5). Small beads (R <C l p ) will 
exclude filaments to a depth of the order of R, and thus 
see an effective compliance larger than that described by 
Eq. (5). Once in this limit, the compliance should again 
scale as This is consistent with our observations: 

the spectra for the smaller beads (0.5 and 0.9 fim) do 
indeed superimpose after scaling (Fig. 8B). 

At low frequencies the slopes of the spectra decrease, 
and the PSDs converge for all beads (Fig. 8A). The onset 
of the slope change, however, does depend on bead size. 
Therefore, for the actin samples it can not be the on- 
set of the elastic plateau as is the case for the PAAm 
data shown below. Instead, we interpret the change 
in slope as the transition, predicted by Eq. (4), from 
a regime of pure shear fluctuations to long-wavelength 
draining modes, for which the spectra should become 
bead-size independent 34 . Such a transition is not seen 
in PAAm because the smaller mesh size inhibits drain- 
ing. The estimate for the crossover frequency for our 
F-actin solutions — on the order of 10 Hz — is also consis- 
tent with the data. 

Fig. 9A-D show the storage and loss shear moduli 
G'(ll>) and G"(u>) obtained for various bead sizes in 
2 mg/ml actin solutions by evaluating the Kramers- 
Kronig integral, Eq. (7). The continuum model should 
apply to the largest beads (Fig. 9D) which are about 20 
times larger than the mesh size. There, scaling is ob- 
served directly in both G' and G" in the range from 
about 5 to 100 Hz. Pure scaling corresponds theo- 
retically (via Eqs. (5) and (7)) to a complex modulus 
G*(uj) oc {iuj) z , so that their ratio is G"/G' = tamtz/2 4 . 
This ratio, averaged between 10 and 100 Hz for the 
largest beads (Fig. 9D), allows a precise estimate of 
z = 0.75 ± 0.02. 

Eqs. (5) and (6) imply scaling of the power spectrum, 
via Eq. (11), as (x^) oc cj -L75 , consistent with the ap- 
parent slope in (Fig. 8A). The scaling exponent z of 
about 3/4 in the macroscopic behavior is distinct from 
the predictions of Rouse-like scaling with exponent 1/2 
observed for flexible polymers (see PAAm below). We 
are unaware of any existing model that predicts z = 3/4 
for the scaling of the complex, macroscopic shear mod- 
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Figure 9: (A-D) Real and imaginary parts of the com- 
plex shear modulus of 2 mg/ml actin solutions, obtained 
from the spectra in Fig. 8 via the Kramers-Kronig rela- 
tion (Eq. (7)) for bead diameters 0.5, 0.9, 2.1 and 5.0 fim. 
Power law scaling of the complex shear modulus emerges 
for the largest beads (D). At low frequencies, an elastic 
plateau is not yet apparent. 

ulus. 

We also note that no plateau is (yet) visible in G'(u>). 
To compare with the literature, we quote values at the 
lowest frequencies we observed. G' at 0.1 Hz was 0.11 ± 
0.03 Pa for 1 mg/ml and 0.17 ± 0.07 Pa for 2 mg/ml 
(actin). Ruddies et a/. 15 measured a value of G' ~ 0.3 Pa 
near the observed onset of the plateau in 0.3 mg/ml 
actin. 

For smaller beads, the apparent G'(u) and G"(u>), cal- 
culated under the assumptions of a single fluid behav- 
ior, show deviations from macroscopic shear elastic be- 
havior (Fig. 9ABC), and a pure scaling regime is not 
reached in the data. This confirms that the assumption 
of an incompressible network is not valid over the whole 
frequency-range of the spectra for the smaller beads. 

Polyacrylamide gels 

Polyacrylamide gels, even at volume fractions close to 
the gelation threshold, are so rigid that we approach the 
limits of our technique. Given the noise floor illustrated 
in Fig. 6, our technique is limited at present to static 
shear moduli of up to a few hundred Pascal. A volume 
fraction of 2% (with relative crosslinker concentration of 
3%) was the lowest concentration that resulted in repro- 
ducibly solidified gels. Fig. 10 shows PSDs of 0.9 fim 



10 



io J 



10' 

' — \ 
N 

I lo- 1 
1 

O 

ft 10 



10" 




Filter 



2.0% 

2.5% 
3.0% 
5.0% 



10" 1 10° io 1 io 2 io 3 io 4 



Frequency (Hz) 

Figure 10: Concentration dependence of PSDs for 
0.9 fxm beads in PAAm gels. The PSD from the 5% gel 
is substantially affected by instrumental noise. The line 
indicates a power law slope of —1.5. The onset of the 
elastic plateau is visible clS cl concentration-dependent 
change in slope at low frequency. 



beads embedded in 2, 2.5, 3, and 5% (volume fraction) 
PAAm gels, linearly averaged over several beads respec- 
tively, and log-binned as described in Materials and Ex- 
perimental Methods. Absence of long-range diffusion of 
the beads confirms that all gels were crosslinked into 
solids. Since the mesh size in the PAAm gels (about 
5 nm in the 2% gels 48 ) is in all cases much smaller than 
the bead size, a continuum elastic model is expected to 
apply throughout. 

PSDs for beads in the softest gels (2%) show a power 
law slope of about —1.5 near 100 Hz, which is distinctly 
less steep than in F-actin. We again observe the slight 
but unexpected downturn in the spectra above about 
3 kHz. The flattening of the spectra below a few Hertz 
reflects the elastic plateau (compare Fig. 11A). This 
plateau extends higher for stiffer gels. For concentra- 
tions above 2.5%, instrumental noise begins to dominate 
the low-frequency part of the spectra (and approaches 
the fixed-bead spectrum, see Fig. 6). Raising the PAAm 
concentration higher above the gelation threshold leads 
to dramatically increased gel rigidity, making the ob- 
servation of fluctuation signals above the noise difficult. 
After calibrating the spectra with the sensitivity factors 
determined for each bead (see Fig. 5), the scatter in the 
spectra typically remained large, unlike for actin. The 
inset in Fig. 13 shows PSDs at 1 Hz for different bead 



sizes in 2.5% gels, indicating this large scatter. 

The 2 and 2.5% gels are not much above the gelation 
threshold and the gel rigidity there changes very rapidly 
with concentration. Therefore, the observed scatter 
is presumably partly due to local gel inhomogeneities. 
Bead polydispersity, as discussed in Materials and Ex- 
perimental Methods, is a further contribution. As for 
actin, we find no correlation with distance from the sub- 
strate surface. The calculated frequency-dependent stor- 
age moduli for the data in Fig. 10 are shown in Fig. 11 A. 
The transformation Eq. (7) has largely eliminated un- 
physical low-frequency noise in these spectra, allowing 
the plateau modulus G'(u — > 0) to be easily estimated. 
For 2, 2.5, 3, and 5% gels the plateau values of G' are 
approximately 2.0, 24, 100, and 280 Pa. 

The effective scaling of G' with concentration is 1.8 
(for concentrations between 3 and 5%), and steeper for 
lower concentrations. Our 2.5% gels thus lie in a regime 
of very rapidly increasing stiffness indeed (this is consis- 
tent with Ref. 49). For larger plateau moduli the upturn 
shifts, as expected, to higher frequencies. The pure scal- 
ing regime beyond the plateau is not reached for any of 
the samples, although its existence is suggested by the 
slope of —1.5 in the PSDs. The corresponding scaling of 
G'(uj) oc ijj 1 ! 2 is consistent with the Rouse model 4 . 

Fig. 11B shows the results of control experiments per- 
formed with conventional cone-and-plate rheology (cour- 
tesy of M. Osterfield, J. Shah, P. Janmey). Samples 
were prepared following the same recipes as for the mi- 
crorheological experiments. Plateau values agree within 
a factor of 2-3. The data show qualitatively the same 
steep dependence on concentration as the microrheolog- 
ical results. Elastic plateaus are clearly visible for 2.5 
and 3% gels whereas for the 2% gels instrumental limits 
are reached, as seen from the large scatter of the data. 

Fig. 12A shows averaged spectra for different bead 
sizes (between 0.5 and 2.1 /im) in 2.5% PAAm gels; 
in Fig. 12B these power spectra are multiplied by the 
respective bead diameters. The i?~ ^dependence in 
Eq. (5) predicts that the rescaled curves have the same 
shape and magnitude. This is approximately true for the 
central frequencies in the spectra, but the high rigidity 
of the PAAm gels causes non-displacement noise to be 
significant both at high frequencies (electronic detection 
noise) and low frequencies (thermal drifts, laser beam- 
pointing and mode instabilities, etc.). The parts of the 
spectra that are dominated by noise show no scaling. 
Due to the bead-size dependence of the detector sensi- 
tivity, the low-frequency noise extends higher for smaller 
beads. 

The storage moduli G'{ui) for three different bead sizes 
2.5% gels are shown in Fig. 13(inset). In contrast to 
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Figure 11: (A) Concentration dependence of the storage 
modulus in PAAm gels, obtained from the spectra in 
Fig. 10. Moduli were calculated for individual beads and 
linearly averaged. Elastic plateaus are clearly visible and 
their extent is found to increase with gel concentration. 
The vertical range is chosen to match (B). (B) Control 
experiments with a cone-and-plate rheometer (Rheomet- 
rics RFS II Fluids Spectrometer) with PAAm gels pre- 
pared in a manner identical to the samples in (A) . Forced 
oscillatory measurements were made at 1% strain and 
at 25°C (data courtesy of M. Osterfield, J. Shah, and 
P. Janmey). 
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Figure 12: (A) Bead-size dependence of PSDs, for beads 
in PAAm gels of 2.5% concentration. Spectra for 0.5, 
0.9 and 2,1 fim bead diameters are averaged over 6, 7 
and 10 beads respectively. (B) Averaged PSDs from (A) 
multiplied by their respective bead diameters. 



12 



T 




Frequency (Hz) 

Figure 13: Variation in the storage modulus of 2.5% 
PAAm gels, for seven different 0.9 /j,m beads. The sub- 
stantial scatter may be partly due to bead polydisper- 
sity. Inset: Values of G' at 1 Hz for multiple beads of 
three sizes. 

what was found for actin, the elastic plateau is clearly 
visible for all samples, and agrees for different bead sizes 
in the same concentration to within about a factor of 
2. The elastic plateau extends to about the same fre- 
quency for all bead sizes. This is in contrast to the low- 
frequency slope changes in actin, which we interpret as 
finite-mesh-size effects, supported by a consistent bead- 
size dependence. With mesh sizes as small as in PAAm, 
draining effects should not play a role, down to the low- 
est frequencies we measure. However, the trend in the 
data is to show apparently smaller moduli for smaller 
beads. We attribute this artifact to a remaining sensi- 
tivity to noise in the integration procedure used to calcu- 
late G' . The detector sensitivity decreases dramatically 
for smaller beads, thus including more noise in the in- 
tegral, leading to the underestimate of G' . The plateau 
moduli estimated by the value of G' at 1 Hz (Fig. 13) 
are 36 Pa (2.1 /an). 24 Pa (0-9 fim), and 17 Pa (0.5 /im 
bead diameter). 

Discussion 

We have introduced the direct use of the Kramers- 
Kronig dispersion integral Eq. (7) to obtain storage 
and loss moduli from power spectra of thermal mo- 
tions. Other methods have been used to obtain shear 
moduli from the thermal mean-square particle displace- 
ment as a function of time: Mason and Weitz give 



an approximate method based on real-valued Laplace 
transforms 28 ; Mason and Wirtz describe a somewhat 
different approximation 29 ' 50 . A limitation of the former 
method is that a specific functional form was postulated, 
fitted to a numerical Laplace transform of particle dis- 
placement data, and then analytically continued to the 
Fourier domain. Possible systematic errors in such a 
scheme are not known. Such transformation methods 
have previously been shown to be sensitive to noise and 
to behavior at the frequency extremes 51 . In contrast 
with Ref. 28, we have applied the transforms described in 
Eqs. (5) and (8) directly to our measured spectra. Both 
of these transform methods are sensitive to the frequency 
extremes, but we believe that the direct transformations 
we have employed are better controlled. 

In order to demonstrate the dependence of our 
Kramers-Kronig transformation on the limits of the in- 
tegration interval, we show in Fig. 14 the transformation 
of a model power spectrum that decreases as a power of 
frequency, above a corner frequency that we arbitrar- 
ily place at 0.1 Hz. We choose the power law to be 
w -i.75^ s i m ii ar to the power law observed in our fluc- 
tuation power spectra for actin solutions. This model 
spectrum corresponds, theoretically, to a complex mod- 
ulus G* that increases as G* oc uj - 75 above the corner 
frequency, and a plateau in the storage modulus G' be- 
low the corner frequency. We show in Fig. 14 the results 
of a transformation carried out with various values of the 
upper and lower cutoffs. The effect of cutoffs is minimal 
in the central portion of the transforms. 

The shear elastic storage and loss moduli (G' and G") 
we found for actin solutions are consistent with other 
experiments 26 ' 27 in the frequency regime where they can 
be compared. Other reported techniques have been lim- 
ited in frequency to a maximum of about 10 Hz. Dis- 
crepancies persist in the literature, 16 and values up to 
300 Pa have been measured for F-actin at a concentra- 
tion of 2 mg/ml. We did not observe changes in vis- 
coelastic behavior in F-actin over periods of at least 24 
hours after polymerization. This suggests that the local 
properties of the solutions probed by the micron-sized 
beads are not affected by slow network changes. On the 
other hand, macroscopic techniques which have observed 
such changes may not have probed the linear response 
regime. 

Actin microrheology by others 26 could be interpreted 
as consistent with u; 3 / 4 scaling, although the authors 
there suggest w 1 / 2 scaling; the frequencies observed were 
close to the onset of the plateau. Macrorheological 
work 15 ' 19 has also reported G* ~ u 1 / 2 . Here too, the 
frequencies studied may have been in the transition be- 
tween plateau and scaling regime. New results from mul- 
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Figure 14: Frequency-cutoff dependence of the Kramers- 
Kronig transformation on a model spectrum. The model 
PSD (not shown) decreases as cj -1 ' 75 above a corner fre- 
quency at 0.1 Hz. The same Kramers-Kronig transfor- 
mation and complex inversion are performed as on the 
actual data in Fig. 9, resulting in the storage and loss 
moduli G' and G" . The lower cutoffs are at 10~ 2 , 1CT 1 - 5 , 
and 10 _1 Hz, and the upper cutoffs are at 10 4 , 10 4 ' 5 , and 
10 5 Hz. For comparison, the curves without any cutoff 
are also included. Deviations from the correct curve ex- 
tend up by about a factor of 2 from the lower cutoff 
and down by about a factor of 10 from the upper cutoff. 
Our actual data in Fig. 9 have an input frequency range 
similar to the curve cut off at 10~ 15 Hz and 10 4 5 Hz. 



tiple light scattering experiments (D. Weitz, personal 
communication) are consistent with a scaling exponent 
close to 3/4 at even higher frequencies than we are able 
to measure. This scaling behavior in semiflexible actin 
solutions may be a universal property of this type of 
network, but is so far unexplained. In particular, there 
is no reason to expect Rouse or Zimm scaling in these 
systems 4 . Bead dynamics in actin networks consistent 
with the power spectra we observe have also been re- 
ported in the limit of very dilute gels (specifically, for 
R — £) 25 ; the- authors suggest a model based on sin- 
gle filament dynamics. We believe that at the (higher) 
concentrations we used, for which the mesh size is sub- 
stantially smaller than bead diameters, our continuum 
elastic approach is correct. 

For solutions of semiflexible polymers, where effective 
entanglement lengths may be much larger than the mesh 
size, 52 ' 53 crosslinking is expected to have a very strong 
effect on the gel rigidity. In cells, the actin cortex is ex- 
tensively crosslinked while the average filament length is 
short. Therefore, a rigidity relevant for cell biology has 
to be measured with a crosslinked actin gel, potentially 
lessening the discrepancies between reported measure- 
ments. 

In both F-actin solutions and PAAm gels, one might 
expect to be able to resolve local inhomogeneities in 
viscoelastic parameters with increasingly smaller beads. 
We did see considerable scatter in the fluctuation power 
spectra, but this spread did not increase for smaller 
beads, as one might have expected. Therefore the in- 
homogeneities in the systems we studied may be smaller 
than effects of bead polydispersity and other instrumen- 
tal errors. More extensive statistical characterization of 
the probe particles and improvement of the instrument 
will reduce the noise, as well as increase the available 
frequency range. These impovements will allow closer 
examination of sample inhomogeneities. 

Our PAAm gels were very close to the gelation thresh- 
old and therefore expected to be inhomogeneous. For 
gels of 2, 2.5, 3, and 5% we indeed find a very steep in- 
crease in G' (Fig. 15). This qualitative behavior and the 
absolute values are consistent with our control experi- 
ments with cone-and-plate rheology and the data match 
literature values in the 3-5% range 49 . Macroscopic meth- 
ods tend to get inaccurate for low shear moduli on the 
order of 1 Pa, whereas our microscopic method under- 
estimates shear moduli above 100 Pa (stiff gels) due to 
noise in the bead position detection. 

In PAAm, the scaling of G*(u) oc uj 1 / 2 is reflected 
by an apparent power law slope of —1.5 in the spectra, 
consistent with the Rouse model 4 . 

A current limitation of our method, so far unex- 
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plained, is a steepening of the power spectrum that we 
consistently observe above a few kHz. We believe this 
is an artifact, because it occurs even for beads in pure 
water. We have ruled out electronic filtering as a cause. 
This effect translates, in the Kramers-Kronig transfor- 
mation of Eq. (7), to a sharp plunge in the calculated 
storage modulus at high frequency, which is of course 
unphysical. Again, this illustrates the subtle nature of 
the frequency sensitivity of the Kramers-Kronig and re- 
lated transformations. Although the downturn in our 
spectra occurs in the vicinity of 10 kHz, its effects are 
apparent in G'(u>) approximately one decade below this. 
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Figure 15: Storage shear modulus of PAAm as a func- 
tion of gel concentration. Solid circles: Microrheology 
results, values estimated from the plateaus in Fig. 11 A. 
Open triangles: Cone-and-plate rheology (data courtesy 
of M. Osterfield, J. Shah, and P. Janmey), values esti- 
mated from the plateaus in Fig. 11B. Stars: Values taken 
from Ref. [49], obtained from dynamic light scattering, 
with a crosslinker fraction of 0.026, similar to our 3% 
preparation. 



Conclusions 

We have used thermal fluctuations of micron-sized 
probes to measure local properties of polymer solutions, 
by a passive and thus non-invasive method, particularly 
appropriate for labile biopolymers. Laser interferome- 
try provides dynamic measurements of probe displace- 
ments with high resolution and up to frequencies not 
previously explored for F-actin solutions. With probes 
that are large compared to the mesh size of the network, 
macroscopic properties such as the complex shear mod- 
ulus G*(lo) can be estimated from the power spectrum 
of fluctuations by using dispersion relations from linear 
response theory. This analysis is made possible by our 
large detection bandwidth. By sampling micron-sized 
volumes we can in principle map out spatial inhomo- 
geneities in soft polymer systems. 

Probes of sizes on the order of the mesh can be used to 
observe deviations from macroscopic behavior, for exam- 
ple the transition to microscopic filament dynamics. Ob- 
serving this breakdown of continuum elasticity may be of 
fundamental importance for understanding the origins of 
elasticity in polymer systems, and has been the subject 
of very recent theoretical attention 54 ^ 56 . Furthermore, 
we have shown that our technique (like other microrhe- 
ological methods using probe particles in polymer solu- 
tions) is predicted to be limited at low frequencies — at 
least in measuring macroscopic shear moduli. This is 
because of expected draining dynamics of network-plus- 
polymer, which is not apparent in macroscopic rheo- 
metric methods. On the other hand, inertial effects of 
both probe and solvent, typically a problem in sensitive 
macroscopic instruments, 1 occur only at MHz frequen- 
cies in microscopic experiments. 

On the micron length scale, a bead in a network should 
behave somewhat like an organelle suspended in the cy- 
toskeleton of a biological cell, and we expect that the 
method will be transferable to measurements in living 
cells, where biological processes are expected to strongly 
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modify viscoelastic parameters. Our technique is cur- 
rently limited by broadband (mainly acoustic) noise to 
rather soft materials, up to shear moduli of a few hun- 
dred Pascal. This limitation will be extended by tech- 
nical improvements. Meanwhile, much remains to be 
done in developing a theoretical understanding of the 
semiflexible polymer systems. 
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